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We introduce and analyze a class of growing geometric random graphs that are invariant under 
rescaling of space and time. Directed connections between nodes are drawn according to influence 
zones that depend on node position in space and time, mimicking the heterogeneity and increased 
specialization found in growing networks. Through calculations and numerical simulations we ex¬ 
plore the consequences of scale-invariance for geometric graphs generated this way. Our analysis 
reveals a dichotomy between scale-free and Poisson distributions of in- and out-degree, the exis¬ 
tence of a random number of hub nodes, high clustering, and unusual percolation behaviour. These 
properties are similar to those of empirically observed web graphs. 


PACS numbers: 89.75.He, 89.75.-k, 64.60.ah 

I. INTRODUCTION 

The theory of random geometric graphs (RGG) en¬ 
ables research of complex networks via geometry m- 
It is attractive to imagine the nodes of a complex net¬ 
work embedded in space, as this can provide useful intu¬ 
ition about an otherwise complicated and abstract dis¬ 
crete structure. The standard RGG can be formulated 
as follows: starting with an empty graph, new nodes ar¬ 
rive at a rate of one per unit time, they are placed at 
random location in the unit square, and attached to all 
existing nodes within distance d. The process is stopped 
when we reach a graph of size N. The properties of net¬ 
works generated this way have been exhaustively studied, 
with particular emphasis on percolation [MS]. In applica¬ 
tions, some specific networked systems are well-described 
by models of this type {e.g. relations between mobile 
telephone users and cell sites (Mil), and the propensity 
of new nodes to attach only to those that are ‘similar’ 
(in the sense of spatial location) mirrors the homophily 
frequently observed in social systems [ 3 ]. However, these 
model networks exhibit a fixed natural scale (see left pan¬ 
els of Fig. and a level of homogeneity that is quite 
atypical of many real world networks m- 

The now textbook explanation for the extreme node 
inhomogeneity observed in many real networks is as a 
result of growth via a mechanism of preferential attach¬ 
ment m- This idea has recently been generalized to 
a geometrical setting in the ‘popularity and similarity’ 
models of web graph formation m ESj. The conse¬ 
quences of preferential attachment for the network in¬ 
clude a scale-free degree distribution and high clustering 
coefficient, however, these same effects can also be gen¬ 
erated by other mechanisms. Recent work in so-called 
‘network cosmology’ considers the modification of ran¬ 
dom geometric graph attachment rules to take account 
of node ‘visibility’ according to the light cone structure of 
some pseudo-Riemannian manifold (e.g., de Sitter space, 
the standard cosmological model [l4Hl6] 'l. which can re¬ 
sult in a scale-free network m- This approach is mean¬ 
ingful in the case of systems with an underlying causal 



FIG. 1: Top row: an RGG in the domain [0,1]^, with N = 
1500 nodes and fixed connection distance radius d = 0.05, 
shown in full (a) and in close-up (b). Bottom row: a scale- 
invariant RGG in the same domain, with N = 10000 nodes 
and influence zone areas Oi = t~^ 12, shown in full (c) and 
in close-up (d). This choice of network ensemble is a special 
case of the general form given in Eq. Q. 


structure, such as citation networks |18) . 

In this article we reverse the question: rather than ask¬ 
ing what causal mechanisms might give rise to a scale-free 
network, we ask what are the consequences of demanding 
that a random geometric graph ensemble is invariant to 
its spatial and temporal scales? We consider a RGG in 
which nodes are created by a unit intensity Poisson pro¬ 
cess in a three-dimensional Euclidean volume (two spatial 
dimensions, plus a time axis) and connected according to 
a function of birth time and spatial proximity. The novel 
aspect we introduce is to demand that the statistics of 
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the resulting graph ensemble depend only on total size 
and not the spatial or temporal scale (see the right panels 
of Fig. [^. As we show, this simple constraint induces a 
range of emergent properties in the generated networks: 
a dichotomy between scale-free and Poisson distribution 
of in- and out-degree, the existence of a random number 
of hub nodes, high clustering, and unusual percolation 
behaviour. 

It turns out that, in this context, scale-invariance is 
synonymous with a particular form of increasing special¬ 
ization of nodes, i.e. that on average the connection ra¬ 
dius is smaller for newer nodes. The increase of special¬ 
ization in general is a common feature of many real-world 
systems; examples of this phenomenon include the pro¬ 
gression of scientific research, the replacement of gener¬ 
alist species by specialists in ecology, and the evolution 
of the world-wide-web. For this last example, one might 
think of nodes as web pages, with the spatial coordinates 
as describing position in some abstract ‘content space’ 
(much like in the models of [IH[I3])j where newer pages 
typically have a smaller scope. Although this story is far 
from a perfect description of the true mechanisms behind 
the formation of web graphs, we will show that some of 
the important properties of the Stanford webgraph can 
be fit to a simple example of our scale-invariant random 
graphs. 

The paper is organised as follows. In Section |ll] we 
present a general formulation geometric random graphs, 
and show how demanding scale-invariance places tight 
constraints on the model specification. The distributions 
of in- and out-degree for scale-invariant geometric ran¬ 
dom graphs are derived in Section m and Section |TV] 
presents some tentative results concerning clustering and 
percolation in these graphs. Finally, in Section we 
briefly illustrate the agreement between certain proper¬ 
ties of scale-invariant geometric random graphs and real- 
world webgraphs. 


II. MODEL DEFINITION 

We consider the following general attachment rule for 
growing random geometric graphs. Nodes appear as a 
unit intensity Poisson process in time (i.e. the time 
between one node appearing and the next is standard 
exponential random variable, independent of everything 
else in the system). Node i is identified by its location 
in space Xi G and time ti G K"*". Each node is as¬ 
signed a spatial zone of influence (or simply ‘zone ’ here¬ 
after) defined as the circular region with center Xi and 
area Ui = f{xi,ti), where / is a positive and piecewise- 
continuous function. For each ordered pair of nodes i 
and j, a directed edge is drawn j —>■ * if Xj lies within 
the zone of i (see Fig. for an illustration). Networks 
are sampled by restricting the model to a finite region S2; 
for simplicity we choose a cylinder with radius R around 
the origin in space and duration T in time. 

For given R and T the expected number of nodes is 



FIG. 2: Illustration of the node attachment rule. In this 
example the nodes arrived in the time-order i,j, k, I and have 
correspondingly ordered zonal areas. 

EN = ttR^T. The standard RGG is represented in this 
class by the choice f{x,t) = tt<P . Note that varying R 
and T but keeping EA^ fixed will cause networks gener¬ 
ated with the standard RGG to have very different sta¬ 
tistical properties; if i? 3> d then the graph will have 
almost no edges, if i? < d then it will be almost fully 
connected. What happens if we demand that the law of 
the random graph depends only on EN and is insensitive 
to the choice of R and T? 

This demand ensures that / must take a particular 
functional form. Gonsidering the law of a single node, 
we have the requirement that for arbitrary scale factors 
A, /r > 0 

f{Xx,g.t) = g{X)h{^j,)f{x,t), (1) 

for some pair of functions g and h. Using polar spatial 
coordinates x = (r, 9 ), it is straightforward to check that 
the only positive continuous solutions to Q have the 
form f{x,t) = a{9)t~^r~^, where a is a bounded func¬ 
tion and / 3,7 are constants. Moreover, it turns out that 
consideration of the graphs as a whole provides a joint 
scaling relation that the exponents j3 and 7 must satisfy. 

For a fixed expected number of nodes EAl, we require 
that the expected number of edges EM does not depend 
on the spatiotemporal scale. The total number of edges 
in a given graph is M = N{k~)/2, where (fc“) is the 
average in-degree of nodes in that graph. For now we 
are only interested in scaling behaviour in large graphs, 
so we assume that (by the law of large numbers) we may 
conflate the arithmetic average over the nodes of a partic¬ 
ular random realization with an ensemble average. Since 
nodes fall as a unit rate Poisson process, the ensemble 
average degree of node i is simply OiT. Ignoring bound¬ 
ary effects and correlations between vertices, we can thus 
compute 

EMRrii / Tf{x,t)dxdt (xT^-^R-"^EN. (2) 

2 Jn 
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Fixing EN as a constant sets R (x T and then EM 
is independent of scale if and only if 

7 = 2(/3-1). (3) 

We thus arrive at the general functional form for deter¬ 
mining the zonal area of a given node: 

(4) 

Following ([^ , the total number of edges scales as EM ^ 
(E7V)^“^, therefore the networks generated are dense if 
/? < 1 and sparse if /3 > 1. Notice also that the depen¬ 
dence on the radial coordinate r means that zonal area 
expands with distance from the origin if /3 < 1 and con¬ 
tracts when /3 > 1. The ensembles of scale-invariant gen¬ 
erated under this rule are therefore not generally trans- 
lationally invariant in space. 

Translational invariance is, however, achieved by the 
special case /3 = l,a(0) = a, that is f{x,t) = a/t. For 
this choice the integral in (H) does not hold, so we must 
check more carefully. Since the maximum degree of a 
vertex is of course less than the number of nodes in the 
graph we can in this case compute 

EM « - / mm{Tf{x,t),EN} dx dt 
2 Jn 

mrR'^ 

= —-— / iiim{T/t,EN} dt (5) 

2 Jo 

cutRJT o; 

= ^ (l + log(EA^)) « ^ENlog{EN). 

So again we find that for this choice of attachment rule 
the total number of edges in the resulting graph is inde¬ 
pendent of the spatial and temporal scale. Fig. con¬ 
trasts a classic RGG with an example of such a transla- 
tionally and spatiotemporally scale-invariant geometric 
random graph. 


III. DEGREE DISTRIBUTIONS 

Let us fix a node i with co-ordinates {xi,ti), and ex¬ 
amine the distribution of in-degree and out-degree. This 
node is the target of an inbound edge from every other 
node inside the cylindrical region 0“ consisting of points 
whose spatial coordinates lie inside the zone of i. Con¬ 
versely, node i has a directed edge leading from it to every 
node j for which Xi G flj (see Fig. I^for an illustration). 
Writing for the region of possible locations of such 
nodes. 


rij = {{x,t) € ^ : tt\x - Xil"^ < f{xi,ti)\, 

^ J (6) 

flf = G n : tt\x - xj'^ < /(a;,t)| . 

The number of nodes that fall within either of these 
regions is a Poisson random variable with mean equal to 



(a) ? (b) t 

FIG. 3: Illustration of the regions Q~ and Q,f, with /3 = 
1/2, Xi = (1.5, 0) in (a) and /3 = 3/2, Xi = (1, 0) in (b). 


the region volume. Therefore the (ensemble) probability 
that a node with known coordinates (xi,ti) has a given 
in- or out-degree is simply 

P(tf =k) = ™Kn.b‘expj-vol(Si.p} 

The spatial nature of the connection mechanism we 
employ induces correlations between nodes, for the sim¬ 
ple reason that if i —>■ j and j —> Z it is more likely that i 
and I are spatially proximate, and hence i ^ 1. However, 
one might hope to obtain a reasonable approximation to 
the degree distribution of a realized network by simply 
integrating Eq. 0 over the uniform distribution of node 
location. For the in-degree p/", this procedure works well. 
We compute 

<*> 

where Lk{x,t) = klog (Tf{x,t)) — Tf{x,t). Laplace’s 
method offers a considerable simplification. Lk achieves 
its maximum in t at 


Uix) = , (9) 

expanding the exponent around this point, we find 

Lfc(a:,t) « fclog(A:) - fc-^ - 1^ . (10) 

Inserting this approximation into © and integrating we 
find 


^P~k 


k^-k 


k^e 


irR^Tkl 


da; dt 


l^fci/2-si//3fc !) /37rl/2i?27^1-l//3 


drdO 


1 

fcl+l//3 


lyajoy^de 


(EN)i/^-i. 


( 11 ) 


where the second line follows from Gaussian integration 
on a time axis extended to ±oo, and the third from using 
Stirling’s approximation for the factorial. 
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FIG. 4: The thin colored lines show the tail distribution of 
in-degree for scale-invariant random geometric graphs with 
/3 = 1/2 (red) and /3 = 3/2 (bine). In each case a single 
graph of size N = 10^ was generated and the empirical tail 
distribution 1 — CDF(fe) calculated. Black lines indicate the 
theoretical result ¥(k i > fc) ~ derived from the power 

law distribution (|11[). 


This calculation becomes exact in the limit fc —>■ oo, 
meaning that the tail of the in-degree distribution has 
a power-law with exponent — (1 + 1//3). One important 
caveat to note is that formally the integrals as they writ¬ 
ten in (111 may not converge, due to singularities at the 
origin in time (/3 < 1) or space (/3 > 1). Mathemati¬ 
cally, this is simply an artifact of us not taking care with 
zones that may overlap the boundary of the sampling re¬ 
gion n, however, this observation is part of a wider phe¬ 
nomenon with major consequences for the distribution of 
out-degrees. 


Examining 0. it is clear to see that nodes arriving 
very early in the process (< ^ 1) enjoy a considerable ad¬ 
vantage over more typical nodes, having extremely large 
zonal areas. When /? > 1, the same is true of nodes close 
to the center of the spatial domain. In fact, any node 
falling in the region 


r = {(r,6l,t) : > Tri?^} (12) 


has a zone the same size as 14. In a particular network 
realization, the possible presence of such a ‘hub’ node 
would have the effect of increasing the degree of almost 
every other node by one, essentially shifting the entire 
out-degree distribution one place to the right. The num¬ 
ber of such hubs is itself a Poisson random variable, with 
mean 


vol(r) = . (13) 


For nodes outside of the hub region in large networks 
(T ^ 1) we can expect that the volume of ilf is dom¬ 
inated by the large T tail, in which we can approxi¬ 
mate the area segment 'k\x — = f{x, t) by the circle 



FIG. 5: Out degree distribution for the scale invariant geo¬ 
metric random graph with /? = 1, a{d) = 1 and ttR^T = 10®. 
The red circles show the out degree distribution averaged over 
100 realizations, while the black curve shows the result of 
equation ( |15[ ). Blue and green triangles are the empirical de¬ 
gree distributions from two of the 100 graphs sampled — their 
horizontal shift results from the random number of hubs. 


7r|a; — = f{xi,t). The result for ,0 ^ 1 is 

i-T 


vol(f4+ \ r) « a{9i)r^ ^ f t ^ dt 

aToCn.e.) 

-1 


= {ENy-^ 




(^-l)(/3-2) /3-1 


(14) 


where the lower limit Tq (r, 0) of the integral was the time 
boundary of the hub region T along the line with hxed 
r and 6. The out-degree distribution of node i is then 
given by a Poisson with mean vol(n)^ \ T), right-shifted 
by a Poisson number of hubs with mean vol(r). 

For the translationally invariant case /3 = 1 and a{9) = 
a there is again a logarithmic correction. The hub region 
for this model is a simp le c ylinder with volume a, but 
the integral in equation (14) evaluates to vol(f4)^ \ T) = 


alog(EiV/a). Averaging over the number h of possible 
hubs, the ensemble out-degree distribution in this case is 




h—0 


{k-h)\ 


(15) 


= ^(a(l + log(EA^/a)))^-“(i+'°s(EW/“)) . 
k\ 


Fig.0 compares this prediction to numerical results. 


IV. CLUSTERING AND PERCOLATION 

Experimentally, we observe that the generated net¬ 
works are typically highly clustered, primarily due to the 
spatial attachment mechanism. There is also an effect 
due to the scale free distribution of in-degree, which can 
bee seen by the following argument. Nodes with small in¬ 
degree typically have small zonal areas and thus nodes 
connecting to them are likely to be spatially close and 
hence connected to each other. The scale-free degree dis¬ 
tribution of in-degree implies that there are very many 
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FIG. 6: Percolation in the scale invariant RGG. (a) shades 
show the fractional size of the strongly connected compo¬ 
nent containing an average node in graphs generated from 
a fixed collection of = 10^ uniformly distributed nodes, 
with various choices of exponent jS and connection length scale 
a{9) = aN^~^. (b) the /? = 1 slice for N = 10^, 10®, and 10^, 
averaged over S = 10®, 10®, and 10 samples, respectively. 


more nodes of small in-degree than large, so they domi¬ 
nate the mean clustering coefficient. Note that this be¬ 
haviour is quite different from that of other networks with 
scale-free distributions, see e.g. Hg. For the translation- 
ally invariant case (/3 = 1, a{0) = a) we can go further 
than this heuristic argument, deriving a scaling relation 
between node degree and clustering. Suppose node i has 
time-coordinate t and node j links to i. If j has time- 
coordinate s and we make the assumption that the over¬ 
lap of the zones of i and j is approximately minja/s, a/t} 
(this is justified if either of s or t is large, which is com¬ 
mon in networks with T ^ 1), then the probability that 
another node I in the influence zone of i is also in the 
influence zone of j is min{l,t/s}. Integrating over pos¬ 
sible values of s, we find C'(t) « t(l -|- log(T/t))/T. For 
the case jS = 1 ,a(9) = a the mean degree of nodes with 
time-coordinate t is simply Tf{x,t) = aT/t, so together 
we find 


C{k) Ki a{l+ log{k/a))/k . (16) 

Another central question in the study of random net¬ 
works is that of percolation, which has been thoroughly 
explored in the standard RGG models [M5]. The pos¬ 
sible existence of infinite connected clusters is by defini¬ 
tion a property of the system in the thermodynamic limit 
of large numbers of nodes. We have numerically inves¬ 
tigated the percolation transition in the scale invariant 
ensemble, for simplicity considering only the rotationally 
invariant case a{0) = aN^~^. The additional scaling 
with N here is required to keep the mean degree fixed as 
N grows and hence obtain a meaningful thermodynamic 
limit. Writing pi for the fractional size of the strongly 
connected component containing node i, we compute the 
size of the largest component p = max^ pi/N, which in¬ 
creases from N~^ to 1 as a varies in [0,oo). The nu¬ 
merical results shown in Fig. [^present a several puzzles. 
Noting the logarithmic scale on the vertical axis of the 
left panel, we can tentatively suggest that for larger (3 
the model appears to percolate for arbitrarily small a 
in the limit N —>■ oo. The inflection point of the level 


lines appears to be close to (0.7,4), shown as a red dot, 
suggesting a critical point in this region. The behavior 
along the slice (3=1, shown in the right panel, is also un¬ 
expected. Our earlier calculation showed that the mean 
degree is l-l-log(Af) in this case, yet it appears that a log¬ 
arithmic correction to a is not required to achieve data 
collapse for different N. 


V. WEBGRAPH COMPARISON 

The properties of scale invariant geometric random 
graphs derived above are not of purely intellectual in¬ 
terest, but in fact show similarity to a particular class 
of real networks: web graphs. These are directed net¬ 
works whose nodes correspond to web pages, and a di¬ 
rected link connects page i to page j if there exists a 
hyperlink on page i referring to page j. Previous au¬ 
thors have proposed causal mechanisms relating to the 
‘popularity’ and ‘similarity’ of pages [HIS] to explain 
the local connectivity properties of web graphs. Here, we 
show that essentially the same behaviour emerges as a 
consequence of the single simple property of scale invari¬ 
ance. We compare the web graph Web-Stanford [50] to 
our scale invariant model. 

The Stanford webgraph exhibits a heavy-tailed in¬ 
degree distribution with exponent approximately 2; we 
therefore must choose /3 = 1 to fit our model. To fit 
the out-degree distribution, we are able to choose the 
angular dependence of influence zone a{9) as we wish. 
For simplicity we take a piecewise constant functional 
form whereby a(0) takes values {ai,... a^} on domains 
of size {2'nwi,... Following the derivation of 

out-degree distribution from earlier we find a mixture 
Poisson 

Ep+ = ^ //c!, (17) 

i 

where q = «;(! -b log(A^/ai)). A reasonable fit 
to the Stanford webgraph is achieved by the values 
{1.4,0.55,0.17} and weights {0.17,0.36,0.47}. In prin¬ 
ciple we could choose a more complex form in order to 
achieve a closer fit, but we find three-component mix¬ 
ture to be sufficient for illustration purposes. From a 
modelling point of view the interpretation is that pages 
in the webgraph can be loosely grouped into three ‘types’ 
that have different typical zones of influence. In Fig. 
the in- and out-degree distributions of the Stanford we¬ 
bgraph is compared with that of a single randomly gen¬ 
erated scale-invariant network, as well as our predictions 
for the ensemble average. 

A simple test of the descriptive power of a fitted model 
of an empirical network is to as if it succeeds in predict¬ 
ing properties that were not explicitly fitted. In Fig. 
we compare the correlation between degree and cluster¬ 
ing coefficient for nodes in our fitted model and the web¬ 
graph. The agreement is by no means perfect, but we do 
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FIG. 7: In-degree (a) and ont-degree (b) distributions of a 
scale-invariant geometric random graph (blue circles) and the 
Stanford webgraph (orange squares), compared with analyti¬ 
cal predictions of Eq. (Ill and Eq. (151 (black lines). 



FIG. 8: Local clnstering coefficient as a fnnction of node in¬ 
degree in the scale-invariant geometric random graph (blue 
circles) and the Stanford webgraph (orange squares), com¬ 
pared with the theoretical prediction of Eq. (161. Note that 
the data has been binned on the fc-axis to improve visibility 
(shown is the average of C{k) for nodes with degree « k). 


observe in both graphs anti-correlation of approximately 
the same magnitude. 


VI. DISCUSSION 

The purpose of this article has been to initiate the 
study of a new class of geometric random graphs that 
are statistically invariant to changes of spatial and tem¬ 
poral scale. We have shown how this simple constraint 
leads to a graph ensemble that is governed by a single 
scale parameter /3 and a function a{6) that determines 
local connectivity and spatial heterogeneity. Analytical 
arguments allow for the prediction of the degree distri¬ 
butions and degree-clustering correlations, which we have 
shown to provide a reasonable fit to those of the Stanford 
webgraph. 

Our model is built on some of the same general ideas 
as explored in the “cosmological” RGGs of e.g. [T7]. One 
might therefore ask if these models are scale-invariant in 
the sense we use here. Certainly, they have a scale-free 
degree distribution (with /3 = 2), however, the spatiotem- 
poral scaling behaviour is quite different. Consider the 
RGG built on (l-l-l)-dimensional de Sitter space. Fix the 
desired number of nodes N « and scale the times- 

pan by a factor A > 0. Then the density 5a should satisfy 
N « i.e. 5a Ne~^'^. In this case the mean 


degree behaves as {k) = e^s5Ads/N « AT, which 
clearly varies with the scaling transformation. Similarly 
the hyperbolic RGG of [H] is also not scale invariant 
with respect to the mean degree — see equation (13) in 
that work. 

The notion of scale invariance is intimately related 
with that of self-similarity of networks, as studied previ¬ 
ously in [22l[23j, for example. In particular, [23] discusses 
an extremely general notion of self similarity in param¬ 
eterized network ensembles. Our model can indeed be 
cast in this form, raising hopes of a theoretical analysis 
of the percolation transition. Unfortunately, none of the 
three specific types of self-similar networks discussed in 
[23] apply directly to our model, so more work is need to 
make this link explicit. 

The work we have presented has several limitations. 
In the model definition, we restricted ourselves to 
two-dimensional space and a single “time” axis that 
parametrizes zonal area. One possible direction for fu¬ 
ture studies would be to generalize to arbitrary dimen¬ 
sions, or to find a conceptually neat way to dispose of the 
time axis. On the subject of the calculations presented 
here, although their results agree well with simulations, 
most have relied on approximations that we have not 
specifically controlled. It would absolutely be desirable 
to formulate more rigorous statements. Finally, some 
readers may find the comparison to empirical webgraphs 
to be somewhat superficial, and they would be correct. If 
scale-invariant geometric random graphs are to find ap¬ 
plications as, for example, null-models in empirical anal¬ 
ysis then a much more nuanced and careful approach will 
be necessary. 

Several more exciting directions for future research 
have also been opened. The most important of these is 
to develop a precise theory percolation in scale-invariant 
graphs. Our numerical work shown in Fig. suggests a 
critical point at a non-trivial coordinate in (a, /?) space; 
clarification of this behavior would be very desirable. 
There is also the possibility that the study these random 
graphs may offer insight into general aspects of real-world 
networks. We have shown here how, for example, scale- 
free degree distributions can occur as a epiphenomenon 
of scale-invariance, which itself is a common property of 
critical systems. The ubiquity of scale-free networks in 
the real-world might therefore point to underlying gen¬ 
erative mechanisms exhibiting self-organized criticality. 
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